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Abstract. In the frame of a three-layer quasi-geostrophic an- 
alytical model of a /-plane geophysical flow, Lagrangian ad- 
vection being induced by the interaction of a monopole vor- 35 
tex with an isolated topographic feature is addressed. Two 
different cases when the monopole locates either within the 
upper or the middle layer are of our interest. In the bottom 
layer, there is a delta function topographic feature, which 
generates a closed recirculation region in its vicinity due to 4o 
the background flow. This recirculation region extends to 
the middle and upper layers, and it plays the role of a topo- 
graphic vortex. The interaction between the monopole and 
the topographic vortex causes complex, including chaotic, 
advection of fluid particles. We show that the model's pa- 
rameters, namely, the monopole and topographic vortices' 
strengths and initial positions, the layers' depths and densi- 
ties are responsible for the diverse advection patterns. While 
the patterns are rather complicated, however, one can single 
out two major processes, which mostly govern fluid parti- 
cle advection. The first one is the variation in time of the 
system's phase space structure, so that within the closed re- 
gion of the topographic vortex, there appear periodically un- 
closed particle pathways by which the particles leave the to- 
pographic vortex. The second one is chaotic advection that 55 
arises from the nonstationarity of the monopole-topography 
interaction. 

Keywords. Three-layer flow, chaotic advection, monopole- 
topography interaction 



1 Introduction 

Generally speaking, topographic vortices are coherent vor- 
tical structures appearing as closed recirculation regions 
over bottom features in the ocean and atmosphere. To- 



pographic vortices play a fundamental role in mass, salin- 
ity and temperature advection in the ocean. Moreover, 
topographic vortices are known to influence the dynam- 
ics of different coher ent structures, such as unrestrictedly 
moving vortices (e.g., Ivan Geffen and Daviesl.ll999l:lpew 



l3; 



bOOllAn and McDonald'. '20051: ICandon and Marshall [20T 
Zavala Sanson et al., 2012). Su ch topographic vortices 
greatly vary in time and size scales ( Baine s and SmithLll993l: 
iBaineSi. il993.) . In this paper, however, we are only inter- 
ested in meso- and synoptic scale topographic vortices due 
to thes e scales are generall y believed to be prevailing in the 



ocean (IChelton et al.Ll201 11) . 



The present paper deals with Lagrangian regular and ir- 
regular (chaotic) advection being generated by a vortex 
monopole interacting with a topographic vortex. The to- 
pographic vortex under investigation is generated by a reg- 
ular t h ree-laye r /-pla ne background flow (e.g., IPedlosky . 



19871: "KozIovI 1 1995k with a delta functi on bottom ir- 



regul a rity within the lower laver (e.g.. Sokolovskiv et all 



1998; Izrailskvetal.. 2004; Kozlov et al 



2005h. 



Then we 



^ '1983'; 



Correspondence to: E. A. Ryzhov 
(ry zho vea @ poi . dvo.ru) 



embed a monopole singular vortex (e.g., Grvanik 
Grvanik and Tevs. 1989: Grvanik et^al, 2 000; Carton, 2001; 
Reznikl I2OIOI: iReznik and Kizneii I2OIOI) either within the 
upper or middle layer. So, these singularities move like 

5 )assive tracers along regul ar backg round flow stream-lines 
Reznik and Kizneii 1200735 l2010h . although, generating a 
complex either periodically (for a time-independent back- 
ground flow) or quasi-periodically (for a periodically time- 
dependent background flow) velocity field in the vicinity of 
themselves. Our main reason for employing such a three- 
layer m odel (e.g., Sokolovskiy, 1997; Ryz hov and Koshej 
l2011a|) is to study Lagrangian advection being induced by 
the monopole-topography interaction when the monopole 
being located within either the upper or middle layers. 
Such monopole positioning can be considered as the sim- 
plest models f or a surface eddy, and f or an inte rthermo- 
cline lens (e.g., ICarton et"al.i. .2002: .Wang and Dewar. .2003: 
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Filyushkin et al. . 2011 ; Filyushkin and Sokolovskiyi 2011 ). 
respectively. However, we should emphasize the paper to 
deal only with surface Lagrangian advection in both cases. 
The middle-layer singular monopole appears on the upper 
layer as a regular vortex, so this configuration generates dif- 
ferent fluid particle advection scenarios due to no singular- 
ities occur in the monopole-topography interaction velocity 
field. 

We investigate two kinds of the monopole-topography in- 
teraction, the first one is an infinity-time interaction and the 
second one is a finite-time interaction. The infinite-time in- 

125 

teraction means the monopole to move in closed regular tra- 
jectories about the topographic vortex's elliptic point for infi- 
nite time due to the constancy of the background flow. How- 
ever, if the background flow depends periodically on time, 
the dynamics of the monopole becomes more complicated. 
So that now the monopole itself can be captured within the 
topographic vortex from the background flow or, on the con- 
trary, be released from the topographic vortex into the back- 
ground flow, and, consequently, be carried away to the infin- 
ity. 

Thus, the main aim of the present study is to investigate 
Lagrangian advection of fluid particles occurring due to the 
velocity field being generated by the infinite and short-term 
monopole-topography interactions . 



95 2 Model formulation 



Substituting Q into O, one can obtain the detailed poten- 
tial vorticity. 



qi = AV^i + /ci(V^2-V^i) + /, 

(72 = AV^2 + (^21 V^l - (^21 - 



(3) 



J^3 



^22) + A:22V^3)+/, 



f'p2 



f'p2 



HigApi^ ^3 - H3gAp2' ^^^l - H2gAp 



^22 



qs = Ai/js - 
where ki - 

H2gAp2 ' 

Since our study concerns only the cases of the upper- and 
middle-layer monopole propagation, we set the lower-layer 
potential vorticity to be always time-independent. However, 
either the upper or middle layer potential vorticity has one 
time-dependent singular value moving with the monopole 's 
center. Hence, we have two sets of singular perturbations of 
the flow. 



Hrr, 



), Qn=qn^ Q3=Q3^ 



(4) 



where m, n = 1,2, m^n, jirn is the monopole's strength 
and is the position of the monopole's singularity within 
the m-layer, q* is the potential vorticity background value. 



and 



with Xi, yi being 



Cartesian coordinates of a fluid particle within the i -layer. 

Potential vorticity ([B should satisfy the potential vorticity 
conservation law in each layer. 



The simplest way to study a quasi-two-dimensional layered g^q- + J(^-^^-) = 0. 
geophysical flow is by exploiting the potential vorticity defi- 
nition in each layer. For a background three-layer flow under 
the ri gid lid approximation, these definitions read (.Pedlosk y. 
1987k 



(5) 



qi Aipi - 



qs Alps 



- /, 92 = A-i/'2 - 



■ — (C2- 



■Ci)- 



To obtain explicit analytical relations for stream-functions 
one can split relations Q]) by making use of the following 
procedure (e.g.. iGryanik and TevsL I1989D . First, we rewrite 
([3]) in a matrix form. 



(1) 



■^ihix,y)- 



145 A* + A* = B, 



<2) + f, 



where i = 1,2,3 corresponds to the upper, middle, and 
lower layer respectively; qi is the z -layer potential vortic- 
ity, Atbi = §^ — ^ is the two-dimensional relative vor- 

' ^ ox oy 

ticity with stream-function t/ji and two-dimensional velocity 
field Ui^Vi; (i, (2 are the interface heights between the upper 
and middle, and the middle and lower layers, respectively; 
h{x,y) = r^(r)) is the Dirac delta function bottom irregular- 
ity with effective volume r; Hi is the z-layer depth; / is the^^° 
constant Coriolis parameter. According to the pressure con- 
tinuity _coiiditioivflie interface heights can be written in the 
form( Pedloskv[ll987h . 



where ^ 



B 




ki 

-(/C21+/C22) k22 

ks -ks ^ 



(6) 



and 



Second, we diagonalize matrix A 



through a similarity transformation, A = SJS~^. Matrix S, 
whose columns are the eigenvectors of matrix A, and diago- 
nal matrix J, whose main diagonal is consisted of the eigen- 
values of matrix A, have the form. 



si = -r-/ s,2 



/(V^3-V^2)P3 



(2) 








-fe(/32-l), 



,(7) 



{g{p2-pi)) ' {g{p3-p2)) ' 

where pi is the i -layer fluid density; g is the gravitational ac- 
celeration; Api = p2 — Pi, and Ap2 = P3 — P2 are the density 
jumps. 



where 



Q^l = -/^22//^21-Q^2//^2l(-/^21-/^22 + fe(<^2-l)), (8) 

a2 = (A:i+A:3 + A:2i+A:22 + Ao)/(2A:3), 



E. A. Ryzhov and K. V. Koshel: Monopole-topography interaction 



3 



-k22/k21 - p2/k21 {-k21 - k22 + ^^^3 (/^2 - 1)) , 

{ki + A:3 + k2i + A:22 - Ao) / (2^) , 



Pi 

P2 



Ao = Y {ki-k^^k2i^k22) -^{-kik^-k^k2i^kik22)- 

Third, we introduce vector $ such that ^ = 
and obtain from © new expression A$ + J$ = 
S~^B, which, taking into consideration relations 

( (^2-p2 Pi- (^1 (^lp2 - OL2P1 \ 

= ^ P2-I I- Pi P1-P2 , and2oo 
y 1 — ^2 ai-l ^2— ai J 
^ = a2 — ai ^ Pi — P2 (^1^2 — (^2^1, has the follow- 
ing detailed form. 



A^2 



-/ - - [(^2 - p2)qi + {Pi - Q^i)<72+ 

7 

+ {(^ip2-^2pi)q3], 



(9), 



ks (a2-l)<^>2 = 

+ (/3l-/32)g3], 
k3{p2-l)^3 = 

+ {a2-ai)q3]. 



-[{P2-I)qi^{l-Pl)q2^ 

7 



--[{l-a2)qi^{ai-l)q2- 

7 



The last step is to obtain explicitly barotropic mode 
and two baroclinic modes ^2, ^3- As it has been men- 
tioned above, we are interested only in the singular pertur- 
bations of form (O (without losing any generality , we put ^^^ 
the b a ckground flow value to be zero, i.e. = dKozlovl 
I995I: llzrailsky et al.l l2004|) ). Hence, by setting boundary 



conditions =0, and d^i/dr\^^^ = to Laplace 

and Helmholtz equations Q, we obtain two sets of Green's 
function superpositions satisfying system Q for the upper 
(m = 1) and middle (m = 2) layer monopole propagation 
cases. 



7 V ^1 

{aiP2-a2Pi)r 



(10) 



2m 



H3 

/^(-ir(/3„-i)Mi 

7 



Hi 



log (ri) 

Ko[Vk3{c^2-l)r*i)- 



3m 



H3 

7 V Hi 
{a2-ai)T 



ko, {a2-l)ri^ 
ii^o(Vfe(/92-lKi)- 



Hz 



where = ^/a 



.2. 



A:3(/32-l)ri 



and m, n = 1,2, m 7^ n. 

Now, introducing a nonvortical plane boundary source flux 
in the form, —Uy, which does not generate any vorticity and 
is compensated by an analogous drain flux (see a detailed 



substantiation in, e.g., llzrailsky et al.l. |2004|) . where is a 
characteristic velocity, one can formulate the final stream- 
functions of the three-layer model with the monopole moving 
within the m-layer. 



-Uy^^im^ai^2m^Pi^ 



3m 7 



(11) 



where as = Ps = l. So, further we will make use of stream- 
functions Sm with the ([TOl) set of functions either for 
the upper (m = 1) or middle (m = 2) layer monopole propa- 
gation case. 

Now, we can introduce certain dimensionless values, 
which will be used further as parameters governing the dif- 
ferent regimes of Lagrangian advection. Introduce length 

]^ /2 

scale L = {ks{a2 — l)) ; velocity scale U; the Rossby 
number, ^ = jj; ; and an effective volume of the topography 
as r = nhoL'^, where h p , L are the height and r adius of an 



corresponding cylinder (ISokolovskiy et al.L 119981) . Then we 



introduce the following governing parameters, 



X-- 



fr 



sHs' 



/Mr) 



(12) 



which characterize the dimensionless topographic vortex 
strength the dimensionless monopole vortex strength, re- 
spectively. Then, by satisfying the quasi-geostrophic re- 
quirement of ^0{e), we set x = ^- Thus, choos- 
ing the following parameters. Hi = 200 m, H2 = 400 m, 
Hs 3000m, pi = 1026.56 kg/m^, p2 = 1027.84 kg/m^, 
ps = 1028.32 kg/m^, we obtain the characteristic horizontal 
topographic vortex scale, L ~ 1.3 • 10^ m. 



3 Equations of motion 

Now, by making use of the dimensionless parameters and 
the geostrophic relations, one can write the equations of mo- 
tion for the monopole 's center and for a fluid particle, being 
advected by the monopole-topography interaction velocity 
field. The monopole motion in the m-layer is governed by 
the following equations of motion. 



dy 



W + X^Kn(C), (13) 



dx 



y^y*rr 



where m = 1,2, and W = W{t) is the dimensionless back- 
ground flow velocity; 



-l{ail32-a2l3i)l+a„,{l3i-l32)Ki{^) + 
7 V C 



+/?m(«2-ai)W-r^^i 

y (0:2 -1) 
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and = \J {x"^)^ + (^m)^ the monopole position with 
m = 1,2 for the upper- and middle-layer monopole motion, 
respectively. System ([T3l) is an elaborated system to govern 
the dynamics of a fluid particle due to the velocity field being 
generated by an exterior background flow intersecting a delta 
function bottom irregularity. Lagrangian advection being de- 
termined by svsten i ([T3]) has been studied recently in the 
frame of barotropic dSokqlovskiy et"allll998l:llzrailskv et al 



-rr 



2004; Koshel and Prants, 2006), two-layer (Kozl ov et al 
2005: Ryzhov and Koshel .201 lb) and three-la yer baroclinic 
geophysical flows (iRyzhov and Koshel 1201 lal) . In our case, 
however, system ([T3]) governs not a fluid particle's motion, 
but a singular vortex's center motion. So, the upper- and 
middle-layer monopoles themselves move as fluid particles 
due to the topographic vortex velocity field. Fluid particles 
of the monopole-topography interaction system, although, 
undergo the joint influence of both the monopole and topo- 
graphic vortex velocity fields. 

Motion of a fluid particle being influenced by the coop- 
erate monopole-topography velocity field obeys to the rela- 
tions. 



dyi 



(14) 



Vi 



dxi 



/ 2 2 

where r*^ = y (x^ - x"^) + {i/i - y^) is the fluid particle 
position relatively to the monopole 's center position, 

260 Pzm{0=^-^({o^n-pn)^ - 1) i^l (0 + 

7 V ^ 



and m, n = 1,2, m 7^ n. 



(^2-1)' 




Wo 
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Fig. 1 : Azimuthal velocities of the topographic vortex within 
the layers. Curves 1,2,3 correspond to the upper, middle, 
and bottom layer, respectively. The horizontal straight line 
indicates constant background velocity value Wo = 0.27r. 



value is the maximal value of the azimuthal velocity in the 
corresponding layer. So, if one chooses the background 
flow to satisfy this condition, then three different- size Tay- 
lor columns will occur due to the bottom irregularity. These 
three columns may be thought of as a discrete Taylor cone. 
Figure [T] depicts azimuthal velocities Vi depending on dis- 
tance r to the topographic vortex elliptic point. We chose 
Wo = 0.2, X = 0.27r to ensure the mesoscale closed regions 
to exist in the all three layers. The points, where the horizon- 
tal line intersects the azimuthal velocity curves, correspond 
to elliptic and hyperbolic critical points of the vortex. Fig- 
ure [2] demonstrates stream-lines of the resulting topographic 
vortex in the upper-layer. The red curve indicates the sepa- 
ratrix dividing the flow into the vortical region and the ex- 
terior flow. Since we also are interested in the middle-layer 
monopole propagation case, the vortical region of the middle 
layer is indicated by the blue dashed curve. 



4 Monopole motion 

First, we briefly analyze system ([T3]) . that governs the 
monopole 's dynamics. An elaborated study of thi s system 
has been conducted in 



2011ah . If the 



iRyzhov and Koshel! 
background exterior flow is constant (W = Wq), system ([T3l) 
is integrable in the sense of the stream-line-trajectory coin- 295 
cidence (e.g., IZaslavskvLll998h . Due to the bottom topogra- 
phy is singular, any nonzero value of Wo always produces 
a closed Taylor column region called a topographic vortex 
within the lower layer. To have such closed regions within 
the middle and upper layers, however, the background ve-300 
locity should be lower than a critical value. This critical 



5 Fluid particle advection 
5.1 Regular monopole motion 

Now we can analyze fluid particle advection being induced 
by the monopole-topography interaction velocity field. Mo- 
tion of a fluid particle is governed by system ([T4b . where the 
right part of the relations comprises the monopole motion 
solution given by ([T3]) . First, we consider the periodic so- 
lution of ([T3]) . This solution, although cannot be expressed 
in an analytical form, is time-dependent with a period be- 
ing equal to the time of the monopole passing a closed tra- 
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Fig. 2: Topographic vortex stream-lines for the upper layer. 350 
The red curve indicates the separatrix. The dashed blue curve 
corresponds to the middle layer separatrix. 



jectory within the separatrices shown in[2l Hence, system 355 
([T4I) is a dynamical system with one and a half degrees of 
freedom, that permits to occur the fluid particle irregular 
dynamics which is conventionally called chaotic advection 
dlref, 1984; Wiggins, 1992; Aref, 2002). Chaotic advection 
manifests itself through expone ntial divergence of close tra- 36o 
jectories in a fin ite tim e (e.g., |Ochtenberg and Lieberman, 
[l983;Zasla vs kyi ll998k . The easiest way to demonstrate the 
chaotic advection manifestation is by constructing Poincare 
sections of system ([T4l) . Figure [3^ shows a Poincare sec- 
tion as hzi = 0.01, yl (0) = —4, corresponding to frequency 355 
a; = 0.1611 of monopole rotation along an orbit shown in 

fig. El 

That half degree of freedom corresponds to a time- 
dependent perturbation, which concerning system (O is the 
monopole motion term comprising strength him- However, 370 
this monopole strength is not the only parameter greatly 
affecting Lagrangian advection, the initial position of the 
monopole is also of great importance. As the initial posi- 
tion parameter, we choose the positions on the ?/-axis due to 
all the stream-lines shown in fig. [Jlintersect this line. Figures 375 
[3]3,c show the system ([T4l) phase space equivalent structures 
for different monopole initial positions, nevertheless, corre- 
sponding to the same stream-line. Hence, the positions on 
the y-8ixis correspond to all the frequencies of the monopole 
rotation about topography. Thus, further we will address how sso 
the monopole 's strength and initial position parameters affect 
the fluid particle dynamics. 

The Poincare section analysis is a very useful technique 
to estimate which part of fluid particles is involved either 
in regular advection or in chaotic advection, however, this 385 
technique fails to show what happens with fluid particles 



in certain moment of time. So, to address the question, 
how these fluid particles move during the monopole pass- 
ing a revolution about the topography, we calculate the num- 
ber of the flow's critical points that appear at e ach instant 
( Ryzhov and KoshelLl2011bl:lRvzhov et al.L 120121) . The sim- 
ple idea of this classification is that the more critical points of 
their initial set survive or, in other worlds, the less topolog- 
ical changes appear during a monopole revolution the more 
regular system ([T4b is. 

5.2 Diagram of the number of the critical points 

As the monopole moves about topography, the number of 
the flow's regular critical points changes, that results in 
the flow t opology altering its characteristics in time (e.g., 
lAref and Brons, 1998). It should be mentioned, that, in the 
upper layer monopole propagation case, one singular critical 
point corresponding to the monopole's center always exists, 
so, we have excluded it from the consideration. Although, 
in the middle-layer monopole propagation case, no singular 
points occur within the upper-layer velocity field due to the 
singular middle-layer monopole appears as a regular one in 
the upper layer. So, making use of the introduced classifica- 
tion, we present diagrams of the number of the regular criti- 
cal points in the upper-layer monopole propagation case and 
in the middle-layer monopole propagation case, respectively, 
depending on monopole's strength k, and initial position y. 
These diagrams depict by color how many regular critical 
points appear at the beginning of the monopole rotation (ini- 
tial critical points) and at the time the monopole passes a half 
of its rotation period (half-period critical points). Figure [4^, 
and fig. lit) correspond to the upper-layer monopole prop- 
agation case and to the middle-layer monopole propagation 
case, respectively. 

Now, we offer a detailed explanation for the diagrams. 
K,<0 region corresponds to counter-rotation of the monopole 
and topographic vortex, and k,>0 corresponds to co-rotation 
of the monopole and topographic vortex. First, we consider 
fig. |4^ depicting the diagram associated with the upper- 
layer monopole propagation case. Figure [5] shows the flow's 
stream-lines at the initial stage of monopole motion and at 
the half-period stage. The red curves are the monopole tra- 
jectories, and the dashed blue curve corresponds to the un- 
perturbed topographic vortex separatrix. 

Also, as a Lagrangian advection measure, we have cal- 
culated the escaping time ( Kozlov and Koshel, 1999, 2000; 



Izrailsky et al.L 12004)) . which is determined as the time a fluid 



particle needs to be carried away by the exterior flow from 
the unperturbed topographic vortex region. This measure is 
an analogue for the Lyapunov exponent and it shows where 
Lagrangian advection progresses faster or slower. Thus, we 
have uniformly distributed within the separatrix 10^ mark- 
ers, and, then, taken into consideration the time they would 
need to cross the line far enough out of the vortex interaction 
(line X = 5). The escaping time distributions are shown in 
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(a) (0, -4) 



(b) (-0.8, -3.3107) 



(c) (0, -0.6936) 



Fig. 3: Equivalent Poincare sections of system (O for the same values of k = 0.01 and the same values of perturbation 
frequency w but different monopole initial positions (xj, y1). 
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-2.5\ 
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-4M 



-5.54 




-1 -0.5 0.5 ^<^ 



-1 -0.5 0.5K^ 



(a) upper-layer monopole (b) middle-layer monopole 
propagation case propagation case 



Fig. 4: Number of the flow's regular critical points by color. 
Blue — 3 initial and 3 or 5 half-period points; purple — 5 
initial and 5 half-period points, grey — 3 initial and 3 half- 
period points red — 3 initial and 1 half-period points; green 
— 1 initial and 1 half-period points. Yellow — 6 initial and^io 
2 half-period points, orange — 4 initial and 2 half-period 
points, brown — 2 initial and 2 half-period points 



fig. [6l where unity of the time is equal to the corresponding 
period of a monopole revolution. A general feature of all the 
subfigures is the almost circle areas of long-live fluid parti- 
cles. These areas correspond to the monopole region, which 
is very intense due to the singularity. Hence, fluid parti cles 
within these areas move mostly regular (IRyzhov and Ko shel, 
I2OI lb), and, therefore, they do not leave the topographic vor- 
tex region. 

The blue color region corresponds to a strong influence 
of the monopole motion. At the initial stage of monopole 
motion, there are three regular critical points to form a hete- 
roclinic structure (see fig. [5^). The topographic vortex can- 
not be distinctly identified due to no hyperbolic point corre- 
sponds to the unperturbed hyperbolic point. So, this case of 
monopole-topography interaction cannot be considered as a 
perturbation of the topographic vortex. Moreover, this ini- 
tial stream-lin e picture resembles a counter-ro tating dipole 



structure (e.g. JVoropayev et al.L I2QQ1I: iRyzhovl 1201 1) . This 



structure changing in time results in that, at the half-period 
stage, there are also three regular critical points to form two 
homoclinic structures each associated either with the topo- 
graphic or with the monopole vortices (see fig. [5j)). Due to 
that topological alteration, fluid particle advection is very ef- 
fective. Most particles are carried away within 5 monopole 
revolutions (see fig. [6^). 

The purple color region corresponds to a moderate influ- 
ence of the monopole motion. At the initial stage, there are 
five regular critical points to form both one heteroclinic and 



E. A. Ryzhov and K. V. Koshel: Monopole-topography interaction 



7 



415 one homoclinic structures (see fig. \5j^). This homocHnic 470 
structure almost coincides with the unperturbed topographic 
vortex separatrix, that indicates this case can be considered as 
a perturbation of the topographic vortex. At the half-period 
stage, the stream-line picture appears as almost the same as 

420 in the aforementioned one of the blue color region (see fig. 475 
[5]l). Lagrangian advection progresses although less efficient, 
but still very fast and it extends on the whole separatrix re- 
gion. No stagnation zone appears within the region (see fig. 

m. 

425 The grey color region corresponds to the least monopole48o 
influence. This positive region differs topologically from 
those presented. Due to the vortices are co-rotating, the ini- 
tial topological structure appears as a co-rotating dipole en- 
veloped by a common topographic vortex separatrix (see fig. 

430 [5t). So, this structure can be considered as a topographic vor- 485 
tex with a double center. Although this double center greatly 
perturbs the fluid particle dynamics, the structure of the to- 
pographic vortex can be revealed during a whole monopole 
revolution (see fig. [Sf). Due to the existence of two always 

435 unbroken centers, fluid particles in the vicinity of the topo- 490 
graphic vortex center move almost regular, however, the sur- 
rounding fluid is carried away very fast (see fig. [S]:). 

The red color region corresponds to a transitional case of 
the monopole-topography interaction. Initially, the stream- 

440 line picture appears as a co-rotating dipole structure (see 495 
fig- Ek). although, during a monopole revolution, the dipole 
structure breaks, so, the singular monopole absorbs the topo- 
graphic vortex elliptic point and becomes a new center of the 
topographic vortex for a certain time. During this time span, 

445 Lagrangian advection within the topographic vortex with the 500 
new singular center is rather regular (see fig. [5]i). How- 
ever, during a whole monopole revolution, almost all the fluid 
from the topographic vortex is carried away (see fig. 

The green color region corresponds to the capturing of the 

450 monopole to be as a topographically trapped vortex with the 505 
singular monopole 's center playing the role of a new topo- 
graphic vortex center. Both at the initial and half-period 
stages, the stream-line portraits comprise only one regular 
critical point that corresponds to hyperbolic point of the to- 

455 pographic vortex. The initial stream-line portrait is shown in 510 
fig. m, while the half-period stream-line portrait appears as 
almost the same as that shown in fig. [5]. Thus, this case can 
be thought of as the topographic trapping of a monopole vor- 
tex. Lagrangian advection, in this case, differs insignificantly 

460 from the case previously addressed (see fig. [6^). 

Now, we consider the diagram shown in fig. ^ for the 515 
middle-layer monopole propagation case. The main differ- 
ence from the upper-layer monopole propagation case is, in 
this case, no singular points appear within the upper-layer 

465 velocity field, so, the monopole vortex appears as a regular 
one and as the regular topographic vortex also can be bro- 
ken. Hence, a merger of the vortices can appear due to both 520 
vortices being regular. This regularity leads to all the half- 
period stream-line pictures to appear almost the same with 



one elliptic point, which is formed by the merger, and one 
hyperbolic point. Also, the lack of a singular point leads to 
that Lagrangian advection is much more regular in compari- 
son with that considered above. 

The yellow color region corresponds to six initial regu- 
lar critical points (see fig. [7^) and two half-period critical 
points (see fig. [TJ)), that both correspond to the topographic 
vortex with the monopole vortex being disappeared since the 
velocity of surrounding flow is too high for a closed circula- 
tion region to be formed. The position of the corresponding 
middle-layer monopole vortex is marked by the cross. Such 
a half-period stream-line portrait is universal for all the color 
regions shown in fig. [TJ). The corresponding escaping time 
distribution is shown in fig. [8^. There is a big stagnation 
region with mostly regular advection corresponding to the 
lower closed region shown in fig. [7^. 

The orange color region is arranged astride the = line. 
This region corresponds to the existence of three initial criti- 
cal points. The difference between the negative and positive 
orange color region initial stream-line portraits is shown in 
fig. [7]:,d. Both the corresponding half-period stream-line 
portraits, however, appear as almost the same as that shown 
in fig. [TJ). Since initially two vortex structures can be reli- 
ably identified, and at half-period stage all these structures 
merge, the escaping time distribution shows very effective 
and intense advection with no stagnation regions progress- 
ing. Figures [8]:,d depict the escaping time distribution in the 
negative and positive k cases, respectively. 

The brown color region corresponds to the existence of 
two initial and half-period critical points. The middle-layer 
monopole does not induce a closed region within the up- 
per layer. Despite that, the middle-layer monopole does 
greatly perturb fluid particle advection. The corresponding 
stream-line portrait does not change topologically during a 
monopole revolution and it appears as almost the same as 
that shown in fig. [Tt. However, on both sides of the = 
line, the advection efficiency is very different. In the < 
zone, advection is very irregular (see fig. [8]l) due to counter- 
rotation of the middle-layer monopole and the topographic 
vortex. On the other hand, in the K:> zone, advection is 
mostly regular, a big stagnation region appears in region of 
the topographic vortex (see fig. [8^), due to co-rotation of the 
middle-layer monopole and the topographic vortex. 

Further, we study how irregular motion of the monopole 
to influence Lagrangian advection 

5.3 Irregular monopole motion 

In this paragraph, we analyze fluid particle advection be- 
ing induced by a non-periodic perturbation consisting of 
periodic background flow oscillation and non-periodic part 
due to monopole irregular motion within the topographic 
vortex. That irregular monopole motion is due to the 
monopole 's singular center moving as a fluid particle in 
the periodically driven velocity field of the topographic 
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(a) initial stage for 
the blue color region 

(ac=-0.5, 2/=-3.5) 




(e) initial stage for 
the grey color region 
(,^ = 0.1,^= -5.5) 




(b) half-period stage for 
the blue color region 
(At =-0.5,2/ = -3.5) 




(f) half-period stage for 
the grey color region 
(^ = 0.1,2/ = -5.5) 





(c) initial stage for the 
purple color region 
(At =-0.1, 2/ = -3.5) 




(g) initial stage for 
the red color region 
(At = 0.3, ^ = -5) 





(d) half-period stage for 
the purple color region 
(At=-0.1,2/ = -3.5) 




for the red color region 
(,t = 0.3,2/ = -5) 



(i) initial stage for 
the green color region 
(At = 0.5,?/ = -4.5) 



(j) half-period stage for 
the green color region 
(At = 0.5,?/ = -4.5) 



Fig. 5: Stream-lines of the flow in the upper-layer monopole propagation case. Red curve corresponds to the monopole motion 
trajectory. Dashed blue curve is the topographic vortex unperturbed separatrix. 



vortex, which is k nown t o produce the irregular dynam- 
ics (e.^..[Sok olovsk iy et al.l, 1998; Kozlov and Koshel, 2001; 54o 
llzrailskv et al.1. I2004 l2008l: Koshel et_aU 120081) . Strictly 
speaking, if one set the background flow to oscillate peri- 
odically. 



= Wo (1 + /iwcosuwt) , 



(15)^ 



where /j^w, and vw are the magnitude and frequency of the 
530 background flow oscillation, then system ([T3l) becomes a sys- 
tem with one and a half degree of freedom permitting the 
chaotic dynamics to occur. Hence, with such an oscillating 
background flow, the monopole can start moving out of the 550 
topographic vortex, and then it can be trapped temporarily 
535 by the topography. And, on the contrary, if the monopole 
starts moving within the topographic vortex, it now can be 
carried away by the exterior flow. It should be mentioned that 
the background flow oscillation also affects the fluid particle 555 



dynamics, resulting in certain number of particles to leave 
the topographic vortex region. However, in our numerical 
simulation, we chose a very small perturbation magnitude 
(fiw = 0.01), so there are very few such particles. So, by 
making use of such a configuration, we study Lagrangian ad- 
vection being mostly induced by the short-term monopole- 
topography interaction. 

Figure [9] depicts an example of fluid particle advection 
being generated by the short-term interaction, while the 
monopole accomplishes a few revolutions within the topo- 
graphic vortex. Figure [9^ shows the initial configuration of 
red and green markers corresponding to the topographic and 
monopole vortex regions, respectively. The unperturbed to- 
pographic vortex region is uniformly filled in with 10^ red 
markers. Also, 1.5- 10^ green markers are placed to dis- 
tinguish the monopole vortex region. The monopole with 
strength tv = 0.1 starts moving out of the topographic vortex 
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Fig. 6: Escaping time distribution in the upper-layer monopole propagation case. 





(a) initial stage for the (b) half-period stage for 
yellow color region the yellow color region 
= -0.7,2/ = -2.5) = -0.7,2/ = -2.5) 





(c) initial stage for the neg- (d) initial stage for the pos- 
ative orange color region itive orange color region 



(^ = -0.5,2/ = -6) 



(^ = 0.5,2/ = -6) 



Fig. 7: Stream-lines of the flow in the middle-layer monopole 
propagation case. Red curve corresponds to the monopole 
motion trajectory. Dashed blue curve is the topographic vor- 
tex unperturbed separatrix. 



(see fig. (9^) at the position with coordinates x = —2, y =590 
—8.4. Then, the monopole vortex is captured by the topo- 
graphic vortex due to chaotic advection (see fig. Next 
fig. [9]: shows the marker distribution as the monopole has 
560 passed a half of rotational period (the black curve points out 



the trajectory of the monopole 's center). A great deforma- 
tion caused by the monopole is clearly seen. Figure [9ji illus- 
trates the particle distribution after the monopole has made 
three whole revolutions about the topography. Few red mark- 
ers from the initial distribution have stayed within the topo- 
graphic region. Last fig. [9^ depicts the monopole leaving the 
topographic vortex region after four revolutions. 

Figure [Tol also depicts a series of marker scattering pat- 
terns, but for the middle-layer monopole propagation case. 
In this case, the monopole starts moving at the position with 
coordinates x = — 1.18,?/ = — 8 and it appears as a regular 
vortex within the upper layer. It results in that a closed recir- 
culation region corresponding to the monopole ceases to ex- 
ist at the half-period stage. Hence, the green markers mostly 
leave the monopole region (see fig. [Tob). However, when 
the closed recirculation region appears again (see fig. [TOji), 
the monopole captures a great deal of the red markers ini- 
tially associated with the topographic vortex. Thus, during 
the topography capturing, the monopole encloses some red 
markers, then after being carried away from the topographic 
vortex, it advects them to the infinity (see fig. [TOfc). 

Figures [9] and [TOl also clearly shows that the particle ad- 
vection is greatly affected to the number of monopole rev- 
olutions about the topography. The longer the monopole 
revolves about the topography the more effective advection 
is. To estimate that short-term monopole influence, we have 
performed a numerical simulation, in which we calculate 
the number of fluid particles escaping the topographic vor- 
tex with respect to the number of the monopole revolutions. 
Since the monopole motion is irregular, two initially close 
monopole trajectories wind the topography very differently, 
with different revolution numbers. Hence, it is impossible to 
predict how many revolutions complete the monopole start- 
ing at a new initial position. Thus, as initial positions for the 
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(e) AC = 0.2, 

7/ = -4 



Fig. 8: Escaping time distribution in the middle-layer monopole propagation case. 











(a) t = 



(b) t = 30 



(c) t = 90 



(d) t = 240 



(e)t = 315 



Fig. 9: Upper-layer monopole propagation case. Particle scattering at the short-term monopole-topography interaction. Red 
and green markers correspond to the topographic and monopole vortex regions, respectively, the blue dashed curve is the 
unperturbed topographic vortex separatrix, and the black curve points out the trajectory of the monopole 's center. Subfigures 
depict markers distribution at the corresponding instant in time. 



595 monopole, we have chosen two intervals of initial positions 
{x = —2, y e [—8.42;— 8.38]) for the upper-layer monopole 
propagation case, and (x = —1.18, G [—8.02;— 7.98]) for 
the middle-layer monopole propagation case. 

Then we have followed the evolution of all the monopoles6i5 
600 starting at these initial position, calculating revolution num- 
ber N of each of those monopoles, and obtained the ad- 
vection efficiency through expression E = na/ni, where 
is the number of advected out of the topographic vortex 
markers, i.e. the markers that have crossed line x = 5, and 620 
605 = 10^ is the initial marker distribution number. It is also 
worth noting, that although some of these monopoles have 
revolved about the topography equal times, Lagrangian ad- 
vection being generated by these monopoles is mostly equiv- 
alent in each case (see fig. [TT]). Indeed, each point in fig. [TT]625 
610 corresponds to one initial position of the monopole. Thus, 



if different initial positions correspond to equal number of 
monopole revolutions TV, then advection efficiency E is suf- 
ficiently similar. 

Figure [TT] depicts advection efficiency E in the upper- 
layer monopole propagation case (see fig. fTTb.b). and in 
the middle-layer monopole propagation case (see fig. fTTb.d). 
By analyzing these subfigures, one can draw several conclu- 
sions. First, TV = 0.5 corresponds to the case of monopole 
passing very close to the topographic vortex but not being 
captured by it. In this case, although, if monopole is very 
weak (a^ = 0.01), it causes a great deal of fluid particle ad- 
vection. A few monopole revolutions are enough for all the 
particles from the topographic vortex region to be carried 
away. Second, the sign of the monopole self-rotation is not 
the main reason of the advection efficiency, but this efficiency 
is mostly determined by Third, evidently, a singular 
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(a)t = (b)t = 18 (c)t = 96 (d)t=120 (e)t = 315 

Fig. 10: The same as in fig. [9] for the middle-layer monopole propagation case. 



monopole (see fig. fTTb.b) causes much more efficient advec- 
tion than a regular one (see fig. fTTb.d). 

6 Conclusions ees 

630 In the frame of a three-layer geophysical flow model, La- 
grangian advection of fluid particle in the vicinity of a 
monopole vortex interacting with a topographic vortex has 
been addressed. Two cases of the monopole propagation 
have been investigated: the upper-layer monopole propaga- 

635 tion, and middle-layer monopole propagation. Such advec- 
tion has been shown to be determined by two most significant 
processes. First, chaotic advection due to the nonstationar- 
ity of the monopole-topography interaction, and, second, the 
appearance or disappearance of closed recirculation zones in 

640 time. Cooperative influence of these processes causes very 
effective Lagrangian advection. Two controlling parameters, 
namely, the monopole's strength and initial position have^so 
been analyzed, and, on the basis of the number of regular crit- 
ical points assessment, a classification of different regimes of 

645 Lagrangian advection has been presented. 

By adding a nonstationary term to the background flow, 
we have analyzed a short-term monopole-topography inter- 
action. If the monopole passes nearly the topographic vortex, 
it still causes a great deal of particles initially located within 

650 the topographic vortex to be carried away. If the monopole 
is captured by the topographic vortex, then it rotates certain 69o 
times about the topography, and, finally, is carried away by 
the background flow. During this passage, the topographic 
vortex almost completely renews its fluid. 

655 Acknowledgements. The reported study was partially supported by 
RFBR, research projects Nos.: 11 -05-00025- A, 11-01-12057-OFl- 
M, 12-05- 31011. 
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